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Abstract 

\ We study a model for the "GSI anomaly" in which we obtain the time evolution of the popu- 

' lation of parent and daughter particles directly in real time, considering explicitly the quantum 

entanglement between the daughter particle and neutrino mass eigenstates in the two-body decay. 

■ We confirm that the decay rate of the parent particle and the growth rate of the daughter parti- 
, cle do not feature a time modulation from interference of neutrino mass eigenstates. The lack of 
\ interference is a consequence of the orthogonality of the mass eigenstates. This result also follows 

from the density matrix obtained by tracing out the unobserved neutrino states. We confirm this 

■ result by providing a complementary explanation based on Cutkosky rules applied to the Feynman 
diagram that describes the self-energy of the parent particle. 
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I. INTRODUCTION 



Recent experiments at the Experimental Storage Ring (ESR) at GSI in Darmstadt 
have revealed an unexpected time dependent modulation in the population of parent ions 
i40pr58+ i42pjn62+ from the Electron Capture (EC) decays i40pr58+ ^ i40Qg58+^^^ 
i42pj^62+ _^ i42j^j62+ _^ [jl]^ ^ phenomenon that has been dubbed the "GSI anomaly." In 
this experiment, changes of the ions' revolution frequencies are detected by the technique of 
time resolved Schottky mass spectrometry. For a small number of stored ions, every decay 
can be resolved. Thus, a time distribution of EC decays of the parent ions can be measured. 
On top of the experimental decay curve, the GSI experiment observed an unexpected time 
modulation with a period of about T ~ 7s. Such a behavior is summarized in Fig. 3-5 of 

Q. 

A theoretical explanation of this remarkable time dependent modulation of the decay rate 
of the parent ion suggests that it is a consequence of the interference between the neutrino 
mass eigenstates in the final state of the two-body decay [11-01 • The authors in Refs. [21, 0, [sj 
argue that the total amplitude of an EC decay is a coherent sum of contributions from 
difference neutrino mass eigenstates. The decay probability is obtained by squaring the 
total amplitude, thus the interference between neutrino mass eigenstates gives rise to the 
observed modulation feature as a consequence of their mixing and oscillations. 

// indeed periodic modulations of EC-decay rates are a consequence of neutrino mixing, 
these experiments would bring an interesting alternative to long-baseline neutrino experi- 
ments for the determination of neutrino mass differences. 

However, this interpretation has been re-examined and criticized in Refs. on the 

basis that it is not the amplitudes that must be summed coherently but the probabilities, 
corresponding to an incoherent addition of the contributions from the different mass eigen- 
states. This approach does not lead to any modulation as the probabilities for the decay 
channels into the different mass eigenstates do not interfere. A similar conclusion is reached 
in ref. l3] by comparing the GSI experiment to other quantum processes, both within quan- 
tum field theory and with quantum mechanical probabilities. 

The theoretical and experimental importance of understanding whether neutrino mix- 
ing and oscillations could lead to time dependent modulations in two body decays where 
neutrinos are a component of the final state warrants an alternative exploration of these 
questions. 

Rather than focusing on any one of these approaches, either summing amplitudes or 
probabilities, in this article, we analyze the two-body decay process differently, by obtaining 
the time evolution of the population of the parent and daughter particles and considering 
explicitly the entanglement between the daughter particle and neutrino mass eigenstates. 



We apply the method developed in Ref. pjj to an analysis of the GSI anomaly, examining 



whether neutrino mixing and oscillations could be responsible for time dependent modula- 
tions in the two-body decay rate. In this approach, we obtain the kinetic equations for the 
populations of the parent and daughter particles directly in real time without the necessity 
to invoke a coherent sum over amplitudes or a sum over probabilities, thereby bypassing the 
theoretical controversy. 

// the time modulation is a consequence of neutrino mixing and oscillations, then this 
phenomenon is robust and does not depend on the details of the parent and daughter nuclei. 
Therefore a simple model of charged current interactions which incorporates neutrino mixing 
but is stripped off the peripheral complications of nuclear matrix elements should describe 
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the essential physical phenomena. 

Therefore, in Section [Tll we introduce a model for the GSI experiment, which captures 
the relevant physical ingredients while neglecting all unnecessary technical complications. In 
Section Hm we obtain the time evolved state emerging from the two-body decay of the parent 
particle. This is a quantum mechanically entangled state between the daughter particle 
and the neutrino, whose time evolution determines completely the number densities of parent 
and daughter particles, unambiguously yielding the time dependence of their population. 
The result of this study confirms that interference between neutrino mass eigenstates is not 
responsible for any modulation in the parent or daughter population, therefore neutrino 
mixing is not the reason behind the GSI anomaly. Section |IV] summarizes our conclusions 
and comments on more recent experimental results. 



II. A MODEL FOR THE GSI ANOMALY 



The EC decays of heavy hydrogen-like particles are governed by charge current weak 
interactions, as shown in Fig. [U^a). If the observed GSI anomaly is a direct consequence 




FIG. 1: (a) exact interaction of EC decays of a parent particle, (b) approximated interaction of 
EC decays in our model. 



of interference between different neutrino mass eigenstates in the final state as proposed 
in Refs. jl, 0, lEI , the technical complications associated with the details of the interaction 
vertices, e.g., spin dependence of fermionic and gauge fields, are irrelevant. In order to 
simplify our calculation, we introduce a bosonic model that captures the main features of 
these EC decays without any unnecessary complication. The two body decay now can be 
represented by the process shown in Fig. [U^b). Our model is specified by the following 
Lagrangian density 

C = Co[M, D] + £oM + Ant[M, D,u,], (2.1) 

with 

Cq[u] = 1 - ^^M^] , (2.2) 

where is a flavor doublet representing the neutrinos 



* = ,? , (2.3) 
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and M is the mass matrix 

(2.4) 



Here, M and D represent the parent and daughter particles, respectively. Their free La- 
grangian density is specified by Cq[M,D]. Also, we consider the simple case of only two 
neutrino fiavors. 

The interaction Lagrangian is analogous to the charged current interaction of the standard 
model, namely 

Cint{x,t) = gM{x,t)D{x,t)ue{x,t), (2.5) 

where g is the coupling constant proportional to the Fermi constant Gp- We note that only 
electron neutrinos enter the interaction because we are considering EC decays. 
The mass matrix is diagonalized by a unitary transformation 

U-\e)MU{9) = f 7 M ; UiO) = ( V (2.6) 

^ ' ^ ' y 7712 J \^ — sm COS a J ^ ^ 

In terms of the doublet of mass eigenstates, the fiavor doublet can be expressed as 

U{e) (2.7) 



In particular i/g = cos 6 ui + sin 6 V2 ■ 



III. NUMBER DENSITIES OF THE PARENT AND THE DAUGHTER PARTI- 
CLES 

Let us consider an initial parent particle state \M{k)) at time t = 0. For the GSI 
experiment, the parent ions are produced with a center velocity of 71% of the speed of light, 
and with a velocity spread Af/f ~ 5 x 10~^ The evolution of the number density of 
parent (M) and daughter (D) particles is obtained from 



NM{t) = (M(^)|e^^*al,(fc)aAf(fc)e-^^*|M(^)), 



iHt 



M{k)), 



(3.1) 



where 7i£,{Q,t) is the number density of daughter particles with momentum Q. Here, the 
annihilation and creation operators are in the Schroedinger picture. We note that e~*^* = 
e~*^°*f7(t, 0) and that the number operators commute with the free field Hamiltonian. U (t, 0) 
is the time evolution operator in the interaction picture, namely. 



U{t,0) = T 



where T is the time-ordering operator. 

Expanding U{t, 0) perturbatively, we obtain U{t, 0)|M) 
■ ■ ■ , where 



|*D(t)) 



(1) 









[ dPxi 







(3.2) 

M) + |^z)(t))('^ + |^D(t))('^ + 
M(fi,ti)D(xi,ti)z/e(fi,ti))l |M(fc)), (3.3) 
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and 



= -g^ [ dh [ £xi r dt2 / d^X2 
Jo J Jo 



M(fi, ti)D{xi, ti)z/e(xi, ti)) 



M{X2, t2)D{x2, t2)ye{x2, ^2)) \M{k)) 



(3.4) 

with z/g = cos^z/i +sin6'z/2. Since \'^D{t))^^^ creates one daughter particle and the initial 
state has none, it is clear that to lowest order, the number of daughter particles is 



(3.5) 



The calculation of the parent population is slightly more involved. The first order state 
has contributions from Fock states with zero or two parent particles M, however the state 
with two parent particles does not conserve energy and its phase varies very rapidly in time 
and averages out in short time scales of order of the inverse mass of the parent particle. 
Therefore to obtain a non- vanishing contribution to the parent population we must consider 
the second order state fl3.4p . 

To second order, there are several contributions, but the only one that is relevant is the 
process in which the first vertex at (x2, ^2) annihilates the initial M creating the intermediate 
state with one (D, z/g) entangled pair, while the second interaction vertex at (xi,ti) annihi- 
lates this (Z), z/e) pair in the intermediate state and creates the M, which has non- vanishing 
overlap with \M). This process is depicted in Fig. [2] and is recognized as the self-energy of 
the parent particle. 



>• 




M 



D 

FIG. 2: Self-energy of the parent particle M, the neutrino line corresponds to a propagator of a 
mass eigenstate. 



Thus to lowest order in g, 



iVM(t) = l + 2Re[(M|^^(t))(2) 



(3.6) 
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A. Number Density of Daughter Particles: The Entangled Daughter- Neutrino 
State 

As pointed out previously, for the number density of daughter particles, we only need to 
consider the interaction Lagrangian up to the first order, namely, Eq. f l3.3p . Expanding the 
field z/g in terms of the fields that create or annihilate the mass eigenstates ui and 1^2, and 
carrying out a standard quantum field theory calculation, we obtain 



(1) 



sin^ 



3^ 1 ^'^>^\v2,p,D^) 



COS 6* i(£;M_£;D_s^ )t 



sin((Ef -E|>-fi2,p)|)■ 



M 
k 



E 



D 



sm{{E¥-E9-n,,,)^) 



M 
k 



E- 



D 



^i,p)/2 



(3.7) 



in which the daughter particle and the neutrinos are entangled Momentum conser- 

vation, a consequence of translational invariance manifest in the Lagrangian density (12. ip . 
enforces p + ^= 0, where p and g* label the momenta of the neutrinos and the daughter parti- 
cle. Also, E^^ and E? are the energies of the parent and daughter particles with momentum 

k and q, respectively. Qi^^ and ^2,^- are the energies of the neutrino mass eigenstates with 
momentum p. In other words. 



E 



M 



m 



E. 



D 



q'^ + m%, Qi^jj = J p'^ + mj, Q 



'2,p 



p'^ + ml (3.8) 



In order to manifestly study the time evolution of the populations and possible time 
dependent phenomena resulting from the mixing of mass eigenstates, we keep the finite 
time-dependence explicitly. As is familiar from Fermi's Golden rule, taking t to infinity 
results in replacing 



sm{{Ef-E^-Q,^)l) 



{Ef-E^-n,,^)/2 



n6iE¥ - E^ - fi,,^ 



D 



(3.9) 



t^oo 



which leads to the standard S-matrix result with the energy conservation at each vertex. 
Here i = 1,2 stand for the neutrino mass eigenstates. 

It is straightforward to calculate (13. 5 p with the state (13. 7p . We find 



noiQ, t) 



8V E¥e^ 

k Q 



COS 



2 Q sin^ {{Ef -Eg- Q,)^] ^ ,in2 ^ sin^ [(E|^ -Eg- Q^^] 



{Ef-E^-n^Y/^ 



^ (3.10) 
where fii 2 = — <5P + "^1 25 corresponding to the specific momentum Q of the daughter 
particle. 

The result (I3.10p is a consequence of the orthogonality of the Fock states associated with 
the mass eigenstates. 
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The production rate of the daughter particle is given by 



dnD{Q,t) 
dt 



k Q 



cosH 2 sin [{E. 



M 
k 



E^ 
Q 



gi^2^2sin [iE¥-E 



D 



(E- -Eg- n,: 



no 



M 
k 



Eg - n,) 



(3.11) 

The time scale of the GSI experiment is about 10 — 100 seconds, corresponding to an 
energy uncertainty AE h/t c:^ 10~^^ — lO^^^eV. Therefore, we can safely take the long 
time limit f l3.9p . leading to the a constant production rate of daughter particles, and the 
total number of daughter particles produced as a function of time is given by 



nD{t) 



^n^(Q, t) 
Q 



Ti cos^ e + T2 sin^ 9 



t. 



(3.12) 



where 



2ng' 
8E¥ 

k 



(27r)3E^fii,2 



(3.13) 



ri^2 are the partial widths, while cos^ 6 and sin^ are the probabilities (or branching ratios) 
associated with each neutrino mass eigenstate. 

From the rate fl3.11l) . we see that there is no interference between the mass eigenstates. 
This is a consequence of the orthogonality of the Fock states associated with mass eigenstates. 
The parent particle decays through two channels, either \i>i) or |i/2) with probabilities cos^ 6 
and sin^ 9 respectively, without interference between them. Obviously when the masses of 
the neutrino vanish the result reduces to the "standard model" decay rate, since Fi = F2. 

This argument becomes clearer upon considering the process of disentanglement of the 
state fl3.3p . The entangled state is disentangled by the measurement resulting in the "col- 



lapsed" state [11 



\VD{Q,t)) 



8VE¥E^ 

k Q 



sin 6 



2,P 



COS 6* 



i(E¥-E^-ni)ir: , 



■sm{{E¥-El-n2)l) 

-my 



iEf-Eg 
sin {{E¥ 



where P = k — Q. It is straightforward to confirm that 

Q 



(3.14) 



(3.15) 



Because p) and 11^2 p) orthogonal with each other, there is no interference between 
these two mass eigenstates. The result f l3.12p is obtained in the long time limit. 



We can further confirm our previous result of the number density of daughter particles 
from the density matrix. The entangled state |\E'£)(t))*^^^ is produced from the decay of a 
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parent particle, correspondingly the density matrix describing this entangled state is 



9' 



1 fsin^e 



8VE¥ 

k q 



cos^ 6 



Dq, l^2,p) {Dq, Z^2,i 



+ 



sin 29 



-sm[{E¥-E5-n,,,)\ 
{Ef-E^-n,,,)/2 



t\-\2 



Sin [{E¥ - E9 - ^],,,) I) 1 rsin [{E¥ - e9 - Q,,,) |) 



{El^ - E^ - n,,,) /2 



X 



_jSrrP_, 
e * 4TT ' 



+ 4n * 



(^f-^?-M/2 

Dg-, (Dg-, Z/2,J? 



(3.16) 



where = — m^, and Vt = \/p'^ + {ml + m\)/2 is the average energy. The density 
matrix contains both diagonal terms, which describe the time evolution of the populations of 
the neutrino mass eigenstates, and off-diagonal terms, which display their interference . 
In the GSI experiment, these neutrinos are not measured, therefore, to calculate the number 
density of daughter particles, we trace out the neutrino states, namely 



no{t) = Tr, 



(3.17) 



Only diagonal terms contribute to the trace because of the orthogonality of the neutrino mass 
eigenstates. Therefore, the number density of daughter particles has nothing to do with the 
interference between neutrino mass eigenstates, which is manifest in the off diagonal density 
matrix elements (coherence). This is consistent with the arguments in [s^ stating that the 
GSI experiment must be described by an incoherent sum over different neutrino states. Once 
again, the answer is Eq. (13. lip in the long time limit. 



B. Number Density of Parent Particles 

Now, let us consider the number density of the parent particles, which follows the same 
line of argument. It proves more convenient to calculate dNM{k,t)/dt, for which we find^ 



dNM{k,t) 
dt 



AE¥ 

k 



M 
k 



E 



D 



sin^ 9 



2 n sm 



{E¥ -Eg- ^],)t 



i2nr Eg] [Ef-Eg-n,) 



[EV - E? - ft.) 
(3,18) 



^ Effectively, we are obtaining the Boltzmann equation for the parent particle, neglecting the build-up of 
the population. 
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In the long time limit this becomes 



dNM{k,t) 
dt 



Ti cos^ e + T2 sin^ 9 



(3.19) 



Clearly, dNM{k,t)/dt = —dnD{t)/dt as the decay of the parent population results in the 
growth of the daughter population with the same rate. This is a consequence of unitarity 
and we can see this by substituting (13.151) and (13. 6p into the unitarity condition 

1 = (M(fc)|f/t(t,0)[/(t,0)|M(fc)) = l+«(^^(t)|^z3(t))(')+2Re[(M|^^(t))(2))]+0((^3)^... 

(3.20) 

Although we have used plane waves to describe our parent and daughter particles our main 
result, the lack of interference of mass eigenstates in the final state is a direct consequence 
of the orthogonality of the mass eigenstates, and this generalizes straightforwardly to the 
case of wave packets. In particular, in reference 0], the wave-packet aspect of the parent 
and daughter nuclei is emphasized as an important ingredient to allow the neutrino mixing. 

However, it is straightforward to show how the main results generalize to the case of 
wave-packets: a wave-packet is a superposition of plane wave components, namely, 

|^(Xo,Po;x,t = 0)) = j d'pf{x,Xo;p,Po)e'nP) (3-21) 

at an initial time t = 0. Here, Xq, Pq are the center position and momentum of wave packet, 
respectively, while x, p and Ep are the position, momentum and energy, respectively. The 
function f{x,Xo;p,Po) specifies the wave function of the particle. 

In our calculation, we obtain the time evolution of each plane wave component \p) from 
which it follows that, 

n*(t) = j d^p\f{x,Xo;p,Po)f n^ip,t), (3.22) 

where n^{p,t) is the parent or daughter population for plane waves obtained above. The 
distribution function / just weights each plane wave component. As demonstrated by our 
calculation in section UTTt the populations n^,{p,t) do not feature oscillations because of the 
orthogonality of the mass eigenstates. Therefore, interference between different neutrino 
mass eigenstates does not appear either in the wave packet treatment, as shown in (13.221) . 
Again, this is a consequence of the orthogonality of neutrino mass eigenstates. Obviously, 
this result is independent of whether the parent or daughter particles are described by plane 
waves or wave packets in agreement with reference (l3|. 

A complementary pathway to the same conclusion is provided by the interpretation of 
(I3.19P in terms of the Feynmann diagram depicted in Fig. [2l This also manifestly leads to 
the conclusion of lack of interference between mass eigenstates because the decay rate of 
the parent particle is the imaginary part of the self-energy. Since the correct propagating 
degrees of freedom are the neutrino mass eigenstates, the total self energy is the sum of 
self-energies with the neutrino mass eigenstates in the intermediate state. As a result, the 
usual Cutkosky rules indicate that the total decay width is the sum of the partial decay 
widths into the mass eigenstates without interference. The real time calculations of the 
decay and production rates presented above confirm this result directly from the evolution 
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of the parent and daughter populations. Introducing a complete set of intermediate states 
between 0) and U{t, in the unitarity condition fl3.20p it is straightforward to find the 
same result as obtained from Cutkosky rules, directly in real time. 

Thus, we confirm the analysis of Refs. |6|-[9| that there is no interference of mass eigen- 
states and we conclude that the GSI anomaly cannot be explained in terms of the interference 
of mass eigenstates in the decay. 



IV. CONCLUSIONS AND DISCUSSIONS 

In this article, we re-examine the GSI anomaly within a framework that is different 
from and complementary to previous work of various groups 0-[9|- The controversy in 
the literature on the theoretical analysis of the GSI anomaly mainly focuses on whether 
probabilities must be summed incoherently joj-loj or amplitudes must be summed coherently 
[2|-|5|. We offer a completely different alternative to study this phenomenon: we obtain 
directly the time evolution of the population of parent and daughter particles taking into 
account that the quantum state arising from the decay of the parent particle is an entangled 
state of the neutrino mass eigenstates and the daughter particle. Our method bypasses the 
issue of summing amplitudes or probabilities and exhibits directly the time evolution of the 
parent and daughter populations. 

Recognizing that if the time modulation of the parent and daughter populations is a result 
of interference phenomena between neutrino mass eigenstates, hence a fairly robust conse- 
quence is independent of the complexities of the parent and daughter nuclei, we introduce a 
simple bosonic model that captures reliably the relevant charged current interaction process 
for EC and manifestly includes neutrino mixing. This allows us to extract the relevant as- 
pects without the peripheral complications associated with spinors, nuclear wave-functions, 
etc. We generalize the recent work [ll| to analyze the GSI anomaly by studying the evolu- 
tion of the distribution functions of the parent and daughter particles directly in time. Our 
starting point is the time evolution of the daughter-neutrino entangled state produced by 
the decay of the parent particle. This treatment also allows us to study the dynamics of 
the daughter particle from the density matrix upon tracing the unobserved neutrino degrees 
of freedom. We show that both the decay rate of the parent particle and the production 
rate of the daughter particle do not feature oscillations arising from the interference of mass 
eigenstates in the final state. This is a direct consequence of the orthogonality of the mass 
eigenstates. 

Furthermore, we provide an alternative field theoretical explanation in terms of the imag- 
inary part of the self-energy diagram of the parent particle. The propagator of the interme- 
diate neutrino states is that of mass eigenstates, therefore Cutkosky rules immediately lead 
to the conclusion that the decay rate is an incoherent sum of probabilities of decay into each 
different mass eigenstate (channels), complementing and confirming our previous analysis. 
Simple arguments based on superposition clearly show that a wave packet treatment of par- 
ent and daughter particles yields the same result, namely no time modulation since there 
is no interference between mass eigenstates in the final state, again a direct consequence of 
orthogonality of mass eigenstates. 

Thus our work confirms the result of Refs. j6|-0] that if the GSI anomaly is a real effect, 
it cannot be explained from the interference of neutrino mass eigenstates. 

More recently, independent experimental efforts have addressed the GSI anomaly: in 
Ref. the EC decay of ^^^Re is studied and no modulation of the decay rate is observed. 
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However, this experiment is different from the one at GSI not only because of the very 
short-hved initial state, but also more importantly because the daughter particle moves in 
a lattice and is restricted to transfer crystal momentum to phonons. 

Another EC-decay experiment with ^^^Pm and an earlier EC-decay experiment with 
^'^'^Eu re-analyzed by Vetter et. al. (Tsf did not observe the modulation in the rates reported 
by the GSI experiment. 

Therefore, our work supports the conclusion against an explanation of the GSI anomaly as 
a consequence of neutrino mixing in agreement with previous work [6i-i9j . These theoretical 
results, combined with emerging independent experimental evidence seem to suggest that if 
the GSI time modulation anomaly is a real phenomenon, its cause is probably associated with 
the details of the GSI experiment, other mechanisms such as neutrino spin precession in the 

or 



static magnetic field of the storage ring [161] , hyperfine level splitting[l7[, spin rotation [18 



perhaps internal nuclear degrees of freedom jl9l| of the parent particle in such an experiment, 
but not a consequence of neutrino mixing. 
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